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Abstract 

It is demonstrated that the stationary Veselov-Novikov (VN) and the stationary 
modified Veselov-Novikov (mVN) equations describe one and the same class of surfaces 
in projective differential geometry: the so-called isothermally asymptotic surfaces, ex- 
amples of which include arbitrary quadrics and cubics, quartics of Kummer, projective 
transforms of affine spheres and rotation surfaces. The stationary mVN equation arises 
in the Wilczynski approach and plays the role of the projective " Gauss-Codazzi" equa- 
tions, while the stationary VN equation follows from the Lelieuvre representation of 
surfaces in 3-space. This implies an explicit Backlund transformation between the 
stationary VN and mVN equations which is an analog of the Miura transformation 
between their (l-l-l)-dimensional limits. 
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1 Introduction 



The Veselov-Novikov (VN) equation 



Ut = au xxx + /3u yyy - 3a(vu) x - 3f3(wu) y 

W X = Uy (1) 

y ^ x 

(here a, f3 are arbitrary constants) was introduced independently in [Q] for real and 
complex-conjugate x, y respectively and arises from the compatibility conditions of the 
linear problem 

V X y = UV 

(2) 

Vt = OLV XXX + (iVyyy ~ 3(XVV X ~ 3[3WVy. 

The modified Veselov-Novikov (mVN) equation 

p t = ap xxx + Ppyyy - 2aVp x - 2(3Wp y - apV x - (3pW y 

W x = Up 2 )v (3) 

Vy = UP 2 ), 

was introduced in j|] and is associated with the two-dimensional Dirac operator 

t/)l = pip 2 
if) 2 , = pip 1 

iPl = aiP l xxx + (3iP l yyy - 2fiW ^ - 3ap x ^ 2 x - f3Wy^ - 2apViP 
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(4) 



= U1PI XX + filPyyy ~ 3/3j>^ ~ 2oV ^ ~ ttV^ ~ IfoW ^ . 

In the (l+l)-dimensional limit a = (3 = ^, y = x equations (|l|) and (j||) reduce respectively 
to the KdV and mKdV equations 

ut = u xxx - 6uu x , p t = p xxx - 6p 2 p x , 

related by the Miura transformation u = p 2 + p x . 

In this paper we establish a similar link between the stationary VN and the stationary 
mVN equations, which in the case a = —(5 assume the following forms 
Stationary VN: 

U XXX 3(VU) X = Uyyy 2>(W!l)y 

W X = Uy (5) 

y — u x . 

Stationary mVN: 

PXXX - 2Vp x - pV X = Pyyy ~ 2Wp y ~ pWy 

Wx = l( P 2 )y (6) 
Vy = l{p 2 )x- 
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Let us introduce the linear system 

r xx =pr y + §(V -p y )r 

(7) 

r yy = pr x + \{W -p x )r 

the compatibility conditions of which coincide with (H). Backlund transformation between 
equations (||) and (||) is given by the formulae 

u = p 2 - 2(lnr ) xy 

v = lV -2{\nr ) xx (8) 
w = lW -2(ln r ) yy 

where r is an arbitrary solution of (|7|). It can be viewed as the stationary analog of the 
Miura transformation. We emphasize that although the (1+1)- and the stationary limits 
of VN and mVN equations are related by Backlund transformations, there is no explicit 
link between the full (2+l)-dimensional equations. 

The origin of this Backlund transformation is purely differential-geometric: both sys- 
tems (|5|) and @ arise in projective differential geometry and describe the so-called isother- 
mally asymptotic surfaces which have been a subject of interest of differential geometers 
of the first half of the century. A number of intriguing properties and interesting examples 
of isothermally asymptotic surfaces can be found in the classical textbooks on projective 
differential geometry §], g, Q, g, @. 

In sect. 2 we recall the standard approach to projective differential geometry following 
Wilczynski 0] and derive the stationary mVN equations (||) playing the role of projective 
"Gauss-Codazzi" equations for isothermally asymptotic surfaces. In this approach linear 
system (^) specifies the radius- vector r of the surface. Particular examples of isothermally 
asymptotic surfaces and exact solutions of the stationary mVN equation corresponding 
to them are discussed in the end of sect. 2. These include quadrics, cubics, quartics of 
Kummer, projective transforms of surfaces of rotation and affine spheres. 

In sect. 3 we recall the Lelieuvre representation of surfaces in 3-space and derive for 
isothermally asymptotic surfaces the stationary VN equation Thus stationary VN 
and mVN equations are just two different parametrizations of one and the same class 
of isothermally asymptotic surfaces and the corresponding Backlund transformation (|8|) 
follows directly. 

2 Surfaces in projective differential geometry 

Following Q we define a surface M 2 in projective space P 3 in terms of solutions of a linear 
system 

r xx = pr y + ~(V -p y )r 

(9) 

r yy = qr x + \{W - q x )r 
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where p, q, V, W are certain functions of x, y. Cross-differentiating (^) and assuming 
r , r x, r y, r xy to be independent, we arrive at the compatibility conditions 

Pyyy ~2p y W ~ pWy = fex ~ ~ (jV X 

W x = 2qp y +pq y (10) 
Vy = 2pq x + qp x 

(see also ||, p. 120 and fl9| ). For any fixed solution p, q, V, W of (jlC|) the linear system (|||) 
is compatible and possesses exactly four linearly independent solutions r = (r° , r 1 , r 2 , r 3 ) 
which may be regarded as homogeneous coordinates of a surface in projective space. For 
definiteness one can think of a surface M 2 in the ordinary 3-space with the radius- vector 
R = (r 1 /r° , r 2 /r° , r 3 /r°). Choosing any other four solutions f = (f° , f 1 ,r 2 , f 3 ) of the 
same system (P) we see that the corresponding surface M 2 with the radius-vector R = 
(fiyp ^fZ/r® ,f 3 /f ) is a projective transform of M 2 , so that any fixed solution p,q,V,W 



of equations ( |10|) defines a surface M 2 uniquely up to projective equivalence. Moreover, a 
simple calculation yields 

Rxx — pRy Q-Rx (11) 



R yy — qR x + bR y 

(a = — 2(lnr°) x , b = — 2(lnr°) J/ ) which implies that x,y are asymptotic coordinates on 
M 2 . In what follows we assume that our surfaces are hyperbolic and the corresponding 
asymptotic coordinates x, y are real (it is not a problem to reformulate results that follow 
in the convex situation regarding x, y as complex conjugate). Since equations (|l0|) specify a 
surface uniquely up to projective equivalence, they can be viewed as the "Gauss-Codazzi" 
equations in projective geometry. 

The most important invariants in projective differential geometry are the projective 
metric 

2pq dxdy (12) 

and the Darboux cubic form 

pdx 3 + qdy 3 (13) 

which define "generic" surface uniquely up to projective equivalence. Projective metric 
(12) gives rise to the projective area functional 

J J pq dxdy. (14) 



Linear system (11) can be the starting point for developing affine differential geometry 
of surfaces referred to asymptotic coordinates x,y. Indeed, let R = (i? 1 , R 2 , R 3 ) be any 
3 nonconstant solutions of (11) viewed as the radius- vector of a surface M 2 in 3-space. 



Since any other 3 nonconstant solutions R = (R 1 , R 2 , R 3 ) are related to R through an 
affine transformation R = AR + b where A is a constant 3x3 matrix and b is a constant 
3- vector (R = const is always a solution of (|i"T|)), system (O) specifies a surface M 2 
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uniquely up to affine equivalence. The compatibility conditions of system (11) 

(q x +aq + \b 2 - b y ) x = 2qp y + pq y 

{p y + bp+ \a 2 - a x ) y = 2pq x + qp x (15) 

Qjy b X 

manifest the "Gauss-Codazzi" equations in affine geometry. 



Remark. The map 

Let us consider linear system (|ll 



iT| ) from projective to affine geometry can be inverted, 
satisfying the compatibility conditions (|l5|). In view of 
a y = b x there exists a function r° satisfying a = — 2(lnr°) x , b = — 2(lnr°) y . Introducing 
the 4-vector r = (r° , r° R 1 , r° R 2 , r° R 3 ) we arrive at (^) with V and W given by 



1 1 

V = p y + bp + -a 2 - ax, W = q x + aq+ -b 2 - b y . (16) 



Thus the "inverse" map pi| ) — » (||) from affine to projective geometry is just the trans- 
formation from affine invariants (p,q,a,b) to projective invariants (p, q,V,W) . 

Isothermally asymptotic surfaces in projective differential geometry are specified by 
the constraint 

p = q (17) 

in which case linear system (^) assumes the form (|7|) while it's compatibility conditions 
( |lO|) reduce to the stationary mVN equation @. The name "isothermally asymptotic" 
which is due to Fubini || reflects the property that in asymptotic coordinates x, y Dar- 
boux's cubic form (]i~3| ) becomes isothermic: p(dx 3 + dy 3 ). These surfaces are projective 
analogs of isothermic surfaces in conformal geometry characterized by the isothermicity 
of the metric in coordinates of lines of curvature. 50 years ago the class of isothermally 
asymptotic surfaces was probably among the most popular ever discussed in the context 
of projective differential geometry. Isothermally asymptotic surfaces arise as the focal 
surfaces of special W-congruences preserving Darboux's curves (zero curves of the Dar- 
boux cubic form) and are uniquely specified by the requirement that the 3-web formed by 
asymptotic and Darboux's curves is hexagonal. We refer to standard textbooks Q, Q, 
P|, P], |1C] for further discussion. 



Remark 1. Isothermally asymptotic surfaces may be defined by the equation 

ln^) =0 (18) 



q 



xy 



which is equivalent to (|l7|) in view of the form-invariance of system @ under the following 
transformations: 



x* = f(x), y*=g(y), r* = y/V r > P* = P9 I '(/ 'f , Q* = if' Kaf 

V*(f') 2 = V + S(f), W*(g') 2 = W + S(g) 
where S(-) is the usual Schwarzian derivative, that is 



W> = £-5(fV 
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(transformation properties of V and W suggest their interpretation as projective connec- 
tions along x— and y— asymptotic curves, respectively). So (18) is just the invariant form 
of (|l7|). It is important that fixing the normalization p = q we fix coordinates x,y up to 
affine transformations x* = ax + b, y* = ay + c, a,b,c = const. Thus any isothermally 
asymptotic surface is endowed with a canonical affine structure. In these coordinates Dar- 
boux's curves are just straight lines x + y = const. The lines x + ky = const also make 
geometric sense: they can be defined in an invariant way as those curves on M 2 which 
have constant cross-ratio with the asymptotic and Darboux's directions. 

Remark 2. Equations @ are invariant under the discrete symmetry (p,V,W) — ► 
(—p,V,W). Geometrically, this means that the class of isothermally asymptotic surfaces 
is self-dual. Indeed, two surfaces ([?]) corresponding to (p, q, V, W) and (— p, —q, V, W) are 
projective duals of each other. 

For isothermally asymptotic surfaces projective metric assumes the form 2p 2 dxdy with 
the corresponding area functional / / p 2 dxdy which is the conserved quantity of the mVN 
equation (||). 

Let us list some of the most important examples of isothermally asymptotic surfaces 
with the emphasize on solutions of the stationary mVN equation @ corresponding to 
them. 



Quadrics correspond to the trivial solution p = 0, W = W(y), V = V(x). 

Projective transforms of rotation surfaces are specified by p = p(x + y), W = 

V = \p 2 + c where p is an arbitrary function of (x + y) and c is an arbitrary constant. 
For c > these are indeed projective transforms of surfaces z = f(x 2 + y 2 ), while the 
cases c = and c < correspond to projective transforms of surfaces z = f(x 2 + y) and 
z = f(x 2 — y 2 ), respectively. Travelling-wave solutions p(x + cy) of equation @ correspond 
to surfaces, which are invariant under one-parameter groups of projective transformations. 
In the case c / 1 the function p is no longer arbitrary and can be expressed in elliptic 
functions (compare with [jbj). 

Cubic surfaces are specified by the following additional constraints in 

V = -^{hip) mi + -(hxp)l + -py, W = --{lnp) yy + -{hip)l + -p x (19) 
|J, see also p. 131. With these V, W equations @ imply 

Vp J y Vp V Vp J x Vp 

Integration of these equations for p would provide a 4-parameter family of exact solutions 
of equation (||): indeed, up to projective equivalence cubics in P 3 depend on 4 essential 
parameters. 

The Roman surface of Steiner is a rational quartic in P 3 with the equation 
(x 2 + y 2 + z 2 - l) 2 = ((* - l) 2 - 2x 2 )((z + l) 2 - 2y 2 ) 
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owing it's name to Steiner who investigated this surface in Rome in 1844. Besides quadrics 
and ruled cubic surfaces the Roman surface of Steiner is the only surface in P 3 possessing 
infinitely many conic sections through any of it's points. This result was announced several 
times: by Moutard in 1865, Darboux in 1880 and Wilczynski in 1908 (see H, 1909 for 
historical remarks). The interest to the Roman surface of Steiner in projective differential 
geometry is due to the remarkable construction of Darboux, relating with an arbitrary 
surface M 2 in P 3 and an arbitrary point p on M 2 an osculating Roman surface of Steiner 
which has the fourth order of tangency with M 2 at this point. Analytically, the Roman 
surface of Steiner corresponds to the choice 

V = ~(lnp) xx + I(lnp)2 - h -p y , W = - l -{\np) yy + l -{\np) 2 y - \p x , (20) 
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,2 



(lnp) xy = -p 

(JlOfl, p. 149-150) implying upon substitution in (pi) the following equations for p: 



_ _ 4 

Pxx — "%PPy 



Pyy = -3PPX (21) 
Qnp)xy = \p 2 - 



These can be explicitely integrated: 



P 



2 



9 f'g 



4(/ + 5) 2 

where the functions f(x) and g(y) satisfy the ODE's 

(/? = (ao + 01/ + a 2 / 2 ) 2 , (g' f = (a - a x g + a 2 g 2 ) 2 . 
Here aj are arbitrary constants. Under the transformation (p, V, W) — * (—p, V, W) equa- 



tions (20) transform to (|19|). This means, that the dual of the Roman surface of Steiner 
is a cubic, and hence the Roman surface itself is a quartic of class 3 (|10||, p. 150). 

The Roman surface of Steiner belongs to a broader class of isothermally asymptotic 
quartic surfaces known as 

Quartics of Kummer investigated by Kummer as singular surfaces of quadratic 
line complexes. Around 1870 Kummer himself constructed plaster models of his quartics 
which now belong to the Gottingen collection (pictures of these models and the necessary 
explanations can be found in p] ) . Analytically, the quartics of Kummer are specified by 
the conditions 

V = ^(lnp) xx + 2(lnp)2, W = f(\np) yy + 2(lnp) 2 y , 

(22) 

(lnp) xy = \p 2 

([10], p. 231). Substituting these V, W in (||) we arrive at 

\{p\p 2 )y)y) ={\{P 2 (P%) X ) ■ (23) 

P 1 )y \P Jx 
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With 

2 = 9 f'g 
P 4 (/ + </) 

equations (p3[) can be rewritten in the form 



(2 1! 

20(/ + 9 )(^-^)-40((/?-( 9 ^ 



(here we used the identities d x = f'df, d y = g d g ). Applying to (24) operator d G /df 3 dg 3 
we arrive at 

W)_«) =n 

dp dg* 

implying that (/ ) 3 and (</) 3 are polynomials of the 6-th order in / and g, respectively. 
Coefficients of these polynomials are not independent and can be fixed upon substitution 
in (H): 

(f'f = P(f), (g'f = P(-g) (25) 

where P is an arbitrary polynomial of the 6th order. Calculations we have presented here 
follow ||, p. 66-69. Formulae (25) reflect the uniformizability of Kummer's quartics via 
theta functions of genus 2 ||i~2|| . Since equations (|2^ ) are invariant under the transformation 
(p, V, W) — ► (— p, V, W) the class of Kummer's quartics is self-dual. 
Quartics of Kummer constitute a subclass of 

Projectively applicable surfaces which are characterized by a condition 

{\np 2 ) xy = cp 2 (26) 

for some constant c. Quartics of Kummer correspond to c = |. Condition (26) means that 
the projective metric 2p 2 dxdy has constant Gaussian curvature K = — c. To investigate 
equations (0) with the additional constraint (p6[) we introduce the anzatz 



V = {\np) xx + l -{\np)l + A, W = (\np) yy + ±(lnp) 2 y + B, (27) 
which implies upon substitution in @ the following equations for A,B: 



Ay = l{l-C){p 2 ) X) B x = l{l-C){p 2 )y, 

{P 2 ) X A+ P 2 A X = (p 2 )yB+p 2 By. 



(28) 



For c = 1 equations Q28| ) are satisfied if A = B = 0. The corresponding surfaces are 
improper affine spheres; they will be discussed below. Here we consider the case c ^ 1. 
Introducing F by the formulae 

A x = -A{\np 2 ) x +F, B y = -B{\np 2 ) y + F (29) 
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and writing down the compatibility conditions of p9| ) with (p8|)i, (|28|)2 we arrive at the 
equations for F 

(30) 

F y = 2cV + |(l-c)^(p 2 (p 2 ) a; ) a; 
the compatibility conditions of which coincide with (p3|) . Inserting in (p3|) 



F x = 2c^ 2 + |(l-c)4(p 2 (p 2 ) 



? = c (7+9)2 PI) 

we end up with the same /, 5 as in ( |25|) . For any p given by ([H]) the functions V and 
can be recowered from ( |27| ) where A,B,F satisfy the compatible system (|28|), (p9|), (|30|). 
Thus for any such p there exists 3-parameter family of surfaces which have the same metric 
2p 2 dxdy and the same cubic form p{dx 3 + dy 3 ) and are not projectively equivalent. In 
general, two projectively different surfaces in P 3 having the same projective metric ( |l2|) 
and the same cubic form ( |l3| ) in a common asymptotic parametrization x, y are called 
projectively applicable. One can show that only for isothermally asymptotic surfaces 
M 2 satisfying ( |26| ) does there exist 3-parameter family of projectively different surfaces 
which are all projectively applicable to M 2 (the value 3 is the maximal possible). We 
point out that geometry of the surface M 2 depends crucially on the value of constant 
c. Isothermally asymptotic surfaces possessing only one-parameter families of projective 



applicabilities have been discussed in [14]. 



Affine spheres constitute an important subclass of isothermally asymptotic surfaces 
specified by the following reduction in @: 

v = P^_l(P_A\ w = Pvv_h(Py)\ (32) 



p 2 V P J V 2 V P J 

After this ansatz the first equation in @ will be satisfied identically while the last two 
imply the Tzitzeica equation for p: 

(lnp) xv = p 2 + -, c = const. (33) 
P 

The cases c 7^ and c = correspond to proper and improper affine spheres, respectively. 
With V, W given by (|32|) equations (|7]) possess particular solution r° = ^Jp. With R = r/r° 
equations (|TTJ) assume the form 

p -„f? _&7? 

^"yy = P^x ^p"-^y 

which become the familiar equations for the radius-vector of affine spheres after adding 
the compatible equation 

p 
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The fact that Tzitzeica's equation solves the stationary mVN equation is also reflected in 
the following nonlocal representation of the mVN equation (in case a = —(3 = 1): 



)xy P 



In the recent paper [|T(| Konopelchenko and Pinkall introduced integrable dynamics 
of surfaces in 3-space governed by the VN equation. One can show that isothermally 
asymptotic surfaces can be interpreted as the stationary points of this evolution. So it is 
not surprising that they are described by the stationary mVN equation (which, as will be 
demonstrated below, is equivalent to the stationary VN). 

Similar class of surfaces, described by the stationary mVN equation (with different real 
reduction) arises in Lie sphere geometry [18]. These are the so-called diagonally cyclidic 



surfaces which in Lie sphere geometry play a role similar to that of isothermally asymptotic 
surfaces in projective geometry. 

3 The Lelieuvre representation of surfaces in 3-space 

We consider a surface M 2 S E 3 parametrized by asymptotic coordinates x, y with the 



radius- vector R(x,y) satisfying equations (11) 



Rxx — pRy + aRx 
Ryy = QRx + bRy 



with the compatibility conditions (|lq). 

The Lelieuvre representation of a surface M 2 [15] is defined by the formulae 



R x = v x v x ,„ ., 

n - v*l (34) 

tXy V X Vy 

where a 3-vector v (called the affine conormal) satisfies the equation 

v xy = uv (35) 

for certain potential u(x, y). In formulae (p3) " x" denotes vector product in E 3 . Equations 



fl34|) are compatible in view of (35). Our aim is to show that in the case of isothermally 
asymptotic surfaces potential u satisfies the stationary VN equation (||). 
Assuming i>, v x , v y to be independent we introduce the expansions 



c l v x + (?v y + (?v 
c^Vx + c 5 u y + c e v 



where the coefficients c* are uniquely expressible through p, q, a, b. Indeed, differentiating 
the first equation (|34]) with respect to x and the second with respect to y we arrive at 
c 1 = a, c 2 = —p, c 4 = —q, c 5 = b. The remaining coefficients can be obtained from the 
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compatibility conditions of the equations for v: u = a y +pq = b x +pq, c 3 = p y +pb, c 6 = 
q x + qa, so that the final equations for v assume the form 

v xy = uv 

v xx = av x - pv y + (p y + pb)v (36) 

v yy = ~1 v x + bv y + (q x + qa)v 
(u = a y + pq = b x + pq) . The compatibility conditions of equations (|^) coincide with 

Remark. Introducing r° by the formulae a = — 2(lnr°) x , b = — 2(lnr ) y one can 
easily check that v = (r°) 2 (R x x R y ) satisfies (36). 

In the case p = q of isothermally asymptotic surfaces the radius-vector R satisfies the 
equations 

Rxx = pRy + o,R x 
Ryy = pR'j: + bR y 

with the compatibility conditions 



(p x + ap + \b 2 - b y ) x = l(p 2 ) y 
(p y + bp+ \a 2 - a x ) y = ^{p 2 ) x 

Cty b X 

while equations for the conormal v assume the form 

v xy = uv 

v xx = av x - pvy + (p y + pb)v (37) 

v yy = ~P y x + bv y + (p x + pa)v 
(u = a y + p 2 = b x + p 2 ). It is a direct check that in this case equations (|37| ) imply 



V XXX 2>W X Vyyy 3WVy (38) 

Vy — U x 
W X = Uy 

where v = |V + a x , w = + b y and V, W are projective invariants (^6[), Since ( ^8[ ) is 
just the stationary VN linear problem, the stationary VN equation (||) for u,v,w follows 
directly. Thus any solution p, V, W of the stationary mVN equation generates solution 
u, v, w of the stationary VN equation according to the formulae 

u = p 2 + a y = p 2 + b x 
v = IV + a x 
w = ~W + by 
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where a = — 2(lnr°) x , b = — 2(lnr°) J/ and r° is a solution of (]?]). This gives the desired 
Backlund transformation (||). 

Let us apply transformation (||) to affine spheres choosing r° = y/p which is a solution 
of linear system (0). With this r° formulae ((^), (^), (|33| ) give the following expressions 
for u, u , w: 



p 3 p 3 \ p/ 3 p 3 \ p 

implying that 

(lnu) xy =u j i u =o ' W =o 

3 u 3 u 

solve (H) . This is reflected in the nonlocal representation of the VN equation (in the case 
a = -P = l) 

1 „ „ 1 



ut = ( d x ud 1 -$d x - d y ud x l ^d y ) (it 2 (lnw) 



y S/ ' J y"" J x U 2"VJ V" V-"-",^ 1113 

which was obtained in [^]. 

As it was pointed out in [pCf| , VN equation is covariant under the Moutard trans- 
formation which can be specialized to the stationary case as well. Geometrically, this 
construction should give the known Backlund transformation for isothermally asymptotic 
surfaces which is generated by a W-congruence preserving Darboux's curves (see e.g. @). 
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